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ABSTRACT 


APL primitives are formally defined by APL/360 
functions. The description is formal in two senses: 
primitives are completely and exactly defined for all 
cases, and the functions are executable on APL/360 and 
are hence working models. 


The descriptions can be used to compare and evaluate 
APL implementations in two ways: 


1. Implemented primitives should produce the same 
results as the corresponding definitions. 


2. Any implementation should properly execute the 
definitions. 


INTRODUCTION 


This is a description of the primitives of APL. They are 
defined by APL\360 functions which describe them to the 
approximate extent of the implementation of APL\360. The APL\360 


User's Manual is the principal reference, and familiarity with it 
is assumed. 


The description is formal in two senses: primitives are 
defined completely and exactly for all cases; and the functions 
which form the description are executable and hence are working 
models. Our intent was to describe the primitives of API more 
completely and more rigorously than does the APL\360 User's Manual, 
but we did not intend the description to be documentation for any 
specific implementation. Hence, we have tended to ignore 
machine-dependent and system-design considerations such as library 
structures and the mechanics of function definition. 


This description can be used to compare and evaluate 
implementations in two ways: 


1. With the exception of machine and system dependencies, 
implemented primitives should behave like the correspending 
definitions for the same arguments. 


2. Since the functions forming the descriptions are themselves 
executable APL, any implementation should execute them 
properly. 


APL was chosen as the lanquage for the description because it 
allows short and concise yet complete and precise definitions. It 
is deficient primarily in primitives for constructing and 
manipulating arrays with components more complicated than scalars; 
this deficiency makes it impossible to formally define and simulate 
APL indexing. 


2 APL SYNTAX 
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We have chosen not to include a formal definition of function 
definition and statement parsing because such a definition adds 
detail which is really not required for understanding. Such a 
definition would in fact be a complete APL interpreter. 


APL SYNTAX 


The formal syntax of APL is relatively simple. Essentially, 
a program consists of a sequence of statements which can be parsed 
into simple expressions. Expressions and statements can be 
informally defined as follows: 


C is a numeric or character constant, 

N is an undefined name or a variable or a [J], 
V is a variable, 

E is an expression, 

F is a function, 

P is a primitive function (see note l), 

Ss is a primitive scalar dyadic function, 

H is a statement label, 

L is a semi-colon list. 


An expression, £, has one of the forms: 


C 

H 

V 

F niladic 

FE monadic 

PLE] E (see note 1) 
BEFE dyadic 

E PLE] E (see note 1) 
S/E reduction 
S/LE)] EF 

S\E scan 

S\CEI & 


Bows EF outer product 


ES.S E inner product 


ECL] subscripted (see note 2) 
N+E specification 
VLLI+E 


A statement has one of the forms: 


E (see note 3) 
H:E 


+E branch 
#:7E 
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Notes: 


1. See Table 1 for those primitives for which a subscript has 
meaning. 


2. The semi-colon list £ in V[L] and £[Z] is of the form 
BeBe wacells where the number of semicolons must be “i+ppV and 
“1+top£& respectively. Any of the expressions in the list may be 
elided, implying all permissible values for that subscript 
position in ascending sequence. 


3. When an expression £ does not contain a specification or branch 
as the function of least precedence, it is assumed that [+ was 
elided on the left, and that the value, if any, of the 
expression is to be displayed. 


Order of evaluation. The relative precedence of functions in an 
expression is positional rather than attributive: precedence 
increases from left to right. Parentheses may be used to delimit 
expressions, and the arguments of a function must be evaluated 
before the function can be evaluated. The rule governing the order 
of evaluation within a statement is this: The rightmost function 
whose arguments are available (i.e. have been evaluated or require 
no evaluation) is evaluated. Thus, for example, the commutativity 
of + is maintained in expressions of the form 


(A*Ax2)+A, equivalent to A+A+Ax2 
ee = 3 3 2 1 (evaluation order) 
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The arguments and results of APL functions are scalars or 
arrays of scalars; a scalar is a number or a character. 
Characters are distinguished from numbers in that no arithmetic 
functions are defined on them. Most current implementations 
require that all scalars forming an array be of the same type, but 
this restriction is not essential. In this paper, the only test 
which prevents mixing numbers and characters is the test in 
COMMACHECK. However, no mixed numeric and character constants yet 
exist in APL, and some of the functions here make the (usually 
trivial) assumption that an argument is of uniform type. 


A vector is a sequence of scalars, formed either by writing 
a constant vector, by catenation, or by an expression involving a 
vector. The reshape function, p, can be used to reshape its right 
argument into an array which has the dimensions specified by the 
left argument. The elements of the result are filled in principal 
order (right-most subscript changing fastest) with scalars chosen 
in principal order from the right argument. If the right argument 
is exhausted, it is repeated cyclically. 
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Examples: 
4nd 2 93 
ieee ee Set 


2 non 2) 3 #5 6 


gO 
4 5 6 

The monadic function p returns a vector of the dimensions of 

its argument. If B is the result of an expression, then: 
eB is a vector of dimensions of the axes of B, 
ppB is the number of axes of 8, and 
x/pB is the number of elements in B. 

Specification, the APL primitive function denoted by the 
symbol + forms an association between the right argument and the 
left argument (which must be a name). The result of a 
specification expression is the right argument. If the name to 


the left of the + is indexed, then only the elements designated by 
the indices are affected. 


array, and is the only function currently permitted in the left 
argument of specification. Vector Indexing is the process of 
selecting particular components of a vector. If xX is a vector and 
I is a scalar such that ITeipX, then X{I] selects the scalar element 
of X located by (I=1pX)/X. If I is not a scalar then the elements 
of xX{I] are obtained by evaluating the scalar for each element of 


I. In all cases, pxX[{IJ] is equal to ol. 


The index origin is the value of J which selects the first 


the vector X is x[0]; in l-origin indexing, $& is: X04) “EBLE Osis 


the value of the index origin, then all indexing can be expressed 
in origin zero as x[I-Q]. 


In general, the components of any array A can be selected by 
the expression ALL] where LZ is a semi-colon list containing ppA 
list elements (which may be arrays) separated by —“1+po0A4 
semi-colons. Tf Bete: of the form) S3e03Ks.207% then pA[Z] is 
CpZ:):5\ God DC pK discs 


If all of the elements of Z are scalars, and Mm is the vector 
formed by catenating these scalars (i.e. MoT SK. ce) > athen>; tin 
terms of vector indexing, the indicated component of A is obtained 
by (,4)£(pA)1M-0] (origin zero). 


An item of LZ may be elided, i.e., the four sequences [; 3; 
3] and [] are all valid. When the [th position is vacant, then the 
value 1(9A)(I] meaning a vector of all permissible values in this 
position in ascending order is assumed. When any of the elements 
is not a scalar, the result is determined by applying the scalar 
case to all combinations taking one scalar from each list element. 


When an index appears to the left of specification in the 
form ALLI+V, the conformability requirement is 
A/(pV)=(pI),(pd),( pK), «ee 
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Function Execution. A defined function consists of a header and a 
sequence of statements. Statements are numbered sequentially 
beginning with l, and the header is referred to as statement 0. 
Function headers, local variables, and statement labels are 
described in detail in the APL\360 User's Manual. Functions are 
executed beginning with statement 1 and, at the completion of any 
statement, proceeding with the next in succession. This sequence 
may be altered by a branch statement. A branch statement is a 
statement formed by the character + followed by an expression. The 
result of this expression must be a scalar or a vector, ana if 
non-empty, the first element X must be an integer. If XeaN 
(origin 1), where W is the number of statements in the function, 
then the next statement executed will be statement %. Otherwise, 
the function terminates, and execution resumes at the point where 
the function was invoked. If the result of the expression to the 
right of the + is an empty vector, then the sequence of execution 
is not changed, and execution proceeds with the succeeding 
statement. Execution of a function terminates after the last 
statement of a function unless it is a branch. 


External Appearance. The external appearance of APL, which is 
described in detail in the APL\360 User's Manual, will not be 
treated formally. The notation used in the definitions is APL\360. 
APL\360 has many important concepts which we have not treated 
formally, but which nevertheless are important to the utility of 


an implementation. We have ignored the library organization and 
system commands. One very important notion is visual fidelity. 
The appearance of a typewritten line corresponds as closely as 
possible to the internal representation of that line. Typing 


errors are corrected in such a way that the statement of correction 
and the correction are both clear and legible. 


Errors are handled in APL\360 in away which facilitates 
recovery and interaction. All errors are detected during 
execution, even though it is possible to detect some errors 
earlier. When an error occurs, no action other than suspension of 
execution is taken. The user is then free to examine the situation 
at the point where the error occurred, and may correct the error 
and resume execution if he desires. The place and manner in which 
errors are detected are shown in each function, and the action 
following detection consists of printing the type of error and then 
terminating execution. The functions usually show the detection of 
errors by a sequence of tests rather than by single more 


complicated expressions. This was done in order to preserve as 
much information as possible about the errors to ease possible 
changes in messages. For example, it might be desirable to 


differentiate LEFT DOMAIN and RIGHT DOMAIN errors. 


Arithmetic ang Fuzz. One of the principal underlying assumptions 
of APL is that the arithmetic primitives are defined on the entire 
continuous domain of real numbers, and that arithmetic is exact. 
We have tended to follow this assumption in the formal definitions, 
so that, even though the notation is APL\360, division, floor, 
ceiling, and the relationals are assumed to be exact. 
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The notion of fuzziness has proved to be so useful in hiding 
the minor errors caused by finite precision and inexact 
representation, that we have included it in definitions of floor, 
ceiling, and the relationals. The definitions are such that 
mathematical identities are preserved. For example, A<B always 
implies B2A. 


Programming Conventions. Since the primary intent was to 
communicate the definitions of APL primitives, we have tended to 
disregard execution efficiency. The functions are meant to be 
executed, and in order to provide as much information as possible 
about their execution through the APL\360 tracing facilities, we 
have minimized the use of multiple specifications and have avoided 
the use of specification in a branch statement. The use of 
indexing has been restricted to the indexing of vectors by scalars 


or vectors. 


The primitive functions which are currently considered part 
of APL are shown in Table 1. All of these have been represented by 
functions except for: the simple arithmetic functions + - x +, 
monadic and dyadic p, indexing, specification, and branching. The 
arithmetic functions have their usual mathematical definitions; the 
others are treated informally below. 


Series or continued-fraction approximations for the transcendental 
functions in order to present executable functions. These 
approximations are not particularly accurate and are not 
recommended for implementation. In general, approximating 
functions should be designed and tailored for the specific host 
hardware. 


Loops and Tests. In general, the order in which loops are 
executed, as well as the order in which tests are made, is 
unimportant, For example, when a scalar function is extended to 
the scalar elements of an array, it doesn't matter which elements 
are chosen first, as long as all elements of the result are 
correctly calculated and stored. Thus, unless a loop contains an 
expression which clearly depends on the sequence followed by the 
induction variable, no inherent meaning should be ascribed to the 
order. This is also true for some sequences of statements. 


x t+ 
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Types: 


st MONADIC DYADIC 


APL PRIMITIVE FUNCTIONS 


s 
s 

s 

s 

s 

s 

s 

s 

s 

s 

so mo 

s PITIMES s 

s NoT u 

u s 

u s 

u s 

u s 

u s 

u s 

u s 

u S 

u s 

u s FNE 
mp - = 
m RAVEL COMMA 
mo XGEN XOF 
u ~- TAKE 
u DROP 


GRADEUP 
GRADEDOWN 


COMPRESS 
EXPAND 
ROTATE 
TRANSPOSE 
MEMBER 
DECODE 
ENCODE 
DMD 


REVERSE 
MTRANSPOSE 


SESE eae 


scalar function 

mixed function 

index has meaning and is origin dependent 
function is origin dependent 

primitive to this report, no formal defini- 
tion given 

no definition exists 


TABLE 1. APL Primitive Functions 


VARIABLES 


Dyadic Identity 
Function Element 


Times 
Plus 
Divide 
Minus 
Power 
Logarithm 
Maximum 
Minimum 
Residue 


ceo 


oe 


oO 
5 
ADOWDDDDAD 


*-O-—-—r—-@6* 1 
ome) ie) 
Do 3 so 
COOnmD @ 


Apply 

for 
logical 
arguments 
only 


N 
L 
L 
L 
N 
L 
L 
L 
N 
N 
L 
L 


Not equal 
Greater 

Not less 
Less 

Not greater 


AAWNVRK We BEd < 
DAAADA 


TABLE 2. Identity Elements of 
Dyadic Scalar Functions 


VARIABLES 


Global Variables. The following variables are parameters to the 
execution of APL expressions, and can be examined and modified by 
the user. 


B - the last random number generated in ROLL. An integer such 
that (820)aB<P 


N - the number of bits ignored in comparisons. 
Q - the value of the index origin, a scalar 0 or l. 


The following represent values which are usually determined by 
hardware: 


C = the character set. In APL\360, there are 256 distinct 
characters, Approximately 150 of these have associated 
printing graphics. 


NB - the floating-point number base. 


WL - the number of digits in the floating-point fraction. 

P and @ are parameters of the random number generator. In APL\ 360 
they are 1+2*31 and 7*5 respectively. Generally, P is chosen 
to be the largest prime which can be stored in the machine 
accumulator and @ is a primitive root of P. 
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The 


remaining global variables are used by the definitions to save 


and pass information: 


F - a character scalar which contains the APL symbol corresponding 
to the scalar dyadic function represented by the example 
function F. 

Fi - the identity element of the scalar dyadic function represented 
by F. 

G - a character scalar which contains the APL symbol which denotes 
the scalar dyadic function represented by the example function 
G. 

I - a scalar integer which is used by AXIS to pass a _ subscript 
value to an indexable function. O=pl is used to indicate an 
elided index, and therefore each function which uses a 
function index must set [+10 before terminating. 

Local Variables. The following conventions are used for naming 

local variables: 

Z- the function result. 

A - the left argument. 

B - the right argument. 

I - axis of application. 

RVLA = ,A 

RVLB = ,B 

RVLZ = ,2 

CA = x/A 

CB = x/B 

CZ = x/Z 

LA = pA 

LB = pB 

LZ = pZ 

XLA = (pA)C(T] 

XLB = (pB)(T] 

XLZ = (9Z)(T] 

TCA = x/I+pA 

TCB = x/I+pB 

TCZ = x/I+¥pZ 
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SCALAR FUNCTIONS 


The scalar functions are those functions, like +, which are defined 
on scalars and give a scalar result. 
to arrays by the adverb functions IP, OP, R, SCAN, and SD, 


V Z+ABS B 


ASCALAR FUNCTION Z+\|B 
Z+B[ -B 


Z+MPLUS B 
ASCALAR FUNCTION Z<«+B 
Z+0+B 


Z+NEGATE B 
ASCALAR FUNCTION Z+-B 
Z+0-B 


Z*+SIGNUM B 
ASCALAR FUNCTION Z+xB 
Z+(B>0)-B<0 


Z*RECIP B 
ASCALAR FUNCTION 2Z++B 
Z+1+B 


Z+A MIN B 

ASCALAR FUNCTION: Z+ALB 
Z+A 

+0 IF ASB 

Z+B 


Z+A MAX B 

ASCALAR FUNCTION Z+A[B 
Z+A 

+0 IF A2B 

Z+B 


Z+A AND B 

ASCALAR FUNCTION Z+AAB 
"DOMAIN' ERROR~Ae 0 1 
"DOMAIN' ERROR~Be 0 1 
Z+ALB 


Z+A OR B 

ASCALAR FUNCTION Z+AVB 
‘DOMAIN' ERROR~Ae O 1 
"DOMAIN' ERROR~Be O 1 
Z+ATB 


SCALAR FUNCTIONS 


These functions are extended 


SCALAR FUNCTIONS 


Vv 


v 


Z+A NAND B 
ASCALAR FUNCTION Z+A*B 
Z+~AAB 


Z+A NOR B 
ASCALAR FUNCTION Z+ANB 
Z+~AVB 


Z+NOT B 

ASCALAR FUNCTION Z+~B 
"DOMAIN* ERROR~Be 0 1 
Z2+1-B 


Z+PITIMES B 
ASCALAR FUNCTION Z+oB 
Z+Bx3.141592653589793 


Z+SHRIEK B 
ASCALAR FUNCTION Z+!B 
AUSES:FACT 
‘DOMAIN’ ERROR O#TYPE B 


"DOMAIN" ERROR(B<0)AB=LB 


Z+FACT B 


Z+ROLL B 
ASCALAR FUNCTION Z+?B 


AGLOBAL VARIABLES:B QO P @ 


"DOMAIN' ERROR O#TYPE B 


'DOMAIN' ERROR(02B)vB2LB 


B«P|Bx@ 
Z+Q+lBxBtP 


Z+A BC B;C3ANI3;BNI;CNI 
ASCALAR FUNCTION 2<A!B 
AUSES:FACT 

"DOMAIN' ERROR O#7YPE A 
"DOMAIN' ERROR O#TYPE B 
C+B-A 
ANI+(A<0)AA=LA 
BNI<(B<0)AB=LB 
CNI«(C<0)AC=LC 
"DOMAIN' ERROR BNIAA#LA 
Z+0 


+0 IF(BNIAANIACNI )v(~BNI )AANIVCNI 


+GAMMA IF(A#LA)VBeLB 
A+(|A)L|B-A 
Z+1 


GAMMA :Z+( FACT B)+(FACT A)xFACT C 
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V Z<LN B 
C1] ASCALAR FUNCTION Z+@B 
G29 "DOMAIN * ERROR Bso 
[3] B+B-1 ; 
C4] 2*B1+Bi2+Bi3+(4xB)+44+(4xB)+5+(9xB)6 


VY Z<ETO B 
C1] ASCALAR FUNCTION Z+*B 
024 2+121-B+1+B+2-Bi3+B:2-Bi5+Bi2 


VY Z+A LOG B 
CiJ ASCALAR FUNCTION Z+A@B 
[2] Z+(@B)+@A 


Exponentiation of Negative Numbers. When A<0O in A*B, the function 
RAPPROX is used to determine whether #8 can be approximated by a 
rational number of the form P:Q. If it cannot be, then the result 
is -(|A)*B. If it can, and @ is even, then A*B is not defined. If 
@ is odd and P is even, the result is (|A4)*B, and if both P and @ 
are Odd, the result is -(|A)*B. RAPPROX returns a 2 2 matrix. The 
two elements of the last column specify the parity of P and @ 
respectively, 0 meaning even andl denoting odd. FE'R as 
irrational, then the result of RAPPROX is 2 29 1, 


VY Z+A EXP B 

[1] ASCALAR FUNCTION Z+AxB 
[2] AUSES:RAPPROX 

[3] +S IF A20 

C4] P+, 0 1 +RAPPROX B 


[5] ‘DOMAIN ' ERROR 0=1+P 
[6] Z+(" 1%*1+P)x*BxelA 
7] +0 
[8] S:Z+*Bxe@A 
Vv 


V P+RAPPROX X;N;3E3B;T 
C1] AUSES: FUZZ 
C2] Pea ZFie.= 412 


[3] N<+10 
C4] E+FUZZ 
[5] B+X 


[6] I7%:;>0 IF BSE 
C7] +IR IF O02N+N-1 


Ce] T+1:+B 

[3] X«LT 

[10] P+oP 

[11] PCL32]+PL32]4PL31]42|% 
[12] BeT-X 

[13] E+FUZZ+ExTxT 

[14] IT 


(151 TR:P+ 2 2p 4 
Vv 
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FACT is used by SHRIEK and BC, It is similar to the gamma 
function, but differs in that it is defined for negative integers. 
! is used on lines 15 and 19 to calculate the gamma function in the 
domain from 0 to l. 


VY Z2+FACT Bsi3UuyF 
{ea Z+1 
C2] F+B-\LB 
[3] +NEG1 IF B<0 
a | T+1+F 
[5] U<+1+B 


C6] +L 
[7] WEG1:I+1+B 
C8] U+F 


[9] L:+F IF U=I 
C10] 2Z<ZxI 

[11] I+I+1 

(12] +b 

[13] B:+NEG2 IF B<oO 
[14] +0 IF 0O=F 
(G59) -ZeZxlr 
[16] +0 

[17] NEG2:2+:2 
[18] +0 IF 0=F 
C19] Z+(Zx!F)tP 


VY gtA RES Bat sc 
C1] ASCALAR FUNCTION 2+A\|B 
{2] Z+B 
{3] +0 IF(0=A)v0=B 
{uJ I+Ax” 1%*(xA)#xB 
[5] C+\A 
[6] BL:tk IP 0<€-|2 
Eda Z+Z-I 
8] +L, 
[9] E:+0 IF(xA)=xB 
[10] Z+<Z+A 

Vv 


V ZA CIRCLE B 
U 


ita ASCALAR FUNCTION 2+AoB 
eg "DOMAIN' ERROR~Ac 84115 
[3] +(L7,£6,L5,L4,£3,62,L1,10,L1,£2,53,04,£5,L6,L7)(CA+8] 
4] £L7:2<ATANI B 

[5] +0 

[6] L6:2+ACOSH B 

[7] +0 

[8] L5:2+ASIN B 

[9] +0 

[10] L4:Z+CIRCLEY B 

[11] +0 

[12] L3:Z+ARCTAN B 

{13] +0 

[14] L2:Z2+ARCCOS B 

Ci5a +0 

[16] Li:Z*+ARCSIN B 

C471 0 


[18] LO:Z*+CIRCLEO B 
[19] +0 
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[20] Li:2+SIN B 
[21] +0 2 
[22] L2:Z«COS B 
[23] +0 
[24] £3:Z<TAN B 
[25] +0 5 
[26] L4:Z<CIRCLEY B 
C273) “0 
[28] L5:Z+SINH B 
[29] +0 
[30] L6:2+COSH B 
[31] +0 
[32] £7:Z+TANH B 
Vv 


V Z+ATANH B;B2 
1] ASCALAR FUNCTION Z+ 70B 
2] "DOMAIN' ERROR 1<|B 
3] B2+BxB 
J Z+Bt1-B2+3-4xB2+5-9xB2+7-16xB2:9-25xB2:11-36xB2:13 


VY Z+ACOSH B 
C1] ASCALAR FUNCTION Z+ 60B 
[2] 'DOMAIN' ERROR 1>B 
[3] Z+@B+ 4OB 


V Z*ASINH B;X;X2 
{1] ASCALAR FUNCTION Z<+ 5OB 
[2] Z+(@(|B)+40B)x~1*B<0 


VY Z+CIRCLE4 B 
[1] ASCALAR FUNCTION Z+ 40OB 
C2] 'DOMAIN' ERROR 1>\|B 
[3] Z«+(_1+BxB)*0.5 


VY Z+ARCTAN B;X;X2 
1] ASCALAR FUNCTION Z+ 30B 
C2] X+(|B 
[3] X2<XxX 
C4] Z*X14X2=344xX2=54+9xX2+74+16«X229 
Cs] Z+Zx 1x*B<0 


VY Z+ARCCOS B 
(1) ASCALAR FUNCTION Z+ 20B 
[23 ‘DOMAIN' ERROR 1<\|B ; 
C3] Z+(00.5)- 10|B 
C4) +0 IF O<B 
(5) Z<«(01)-Z 


SCALAR FUNCTIONS THAT USE FUZZ ES 


V Z+ARCSIN B;X3;X2 s 
Es) ASCALAR FUNCTION Z+ 10B 
[2] 'DOMAIN' ERROR 1<\|B 
C3 X2<XxX , ; 
C4] Z*((1-X2)*0.5 )xXF1-2xXQE9-2XXQEH-12XK2ET-12xK2z9 


V Z<CIRCLEO B 
C1] ASCALAR FUNCTION Z+00B 
[2] "DOMAIN' ERROR 1<\|B 
[3] Z+(1-BxB)*0.5 


V Z<«SIN B;C 


{1] ASCALAR FUNCTION Z+10B 
[2] C#14+2x 14115 
[3] Ze-/(B«C)+!C 


V 2+COS B 
Ea) ASCALAR FUNCTION 2Z+20B 
[2] Z+10(00.5)-8 


V Z+TAW B 
[ij ASCALAR FUNCTION Z+30B 
[2] Z<«(10B8)+208 


VY Z«CIRCLE4 B 
4] ASCALAR FUNCTION 2Z+4OB 
C2] Z+(1+BxB)*0.5 


VY Z2<SIWH B 
{1] ASCALAR FUNCTION Z2<5OB 
[2] Z+0.5x-/*B,-B 


| 


aj 


Z+COSH B 
[1] ASCALAR FUNCTION Z<60B 
2a Z+0.5x+/*B,-B 


i B 
AR FUNCTION Z+70B 
B)+60B 


rs wi 


t TH + 
>) 
wm 
Ot 


SCALAR FUNCTIONS THAT USE FUZZ 


The primary use of fuzz occurs in floor, ceiling, and the 
relationals. The first two functions below define absolute and 
relative fuzz, respectively. The three global variables used have 
the following meanings: 


48 The base used to represent the floating-point fraction. In 
System/360, the base is 16. 


WZ The number of digits, base WE, forming the floatino-point 
fraction. In System/360, the number of digits in the fraction 
is 14, 


16 


I= 


C1) 
C2] 
C3] 
C4] 


SCALAR FUNCTIONS THAT USE FUZZ 


The number of bits, or binary digits, to be ignored in 
comparisons. As the functions demonstrate, this notion is 
valid on machines which do not use binary encoding, and 


Wilkinson error analysis techniques are valid. 


V Z+FUZZ 
AGLOBAL VARIABLES:N NB WL 
Z+((2*N)-1)*xNB*-WL 

Vv 


V Z+A RFUZZ B 
AGLOBAL VARIABLES:NB 
AUSES:FUZZ 
ARELATIVE FUZZ 
Z+FUZZxNBx[NBe(|A)T IB 

Vv 


True ceiling and true floor have been included for illustration and 
emphasis in the other functions. 


VY Z+TCL B 
ATRUE CEILING 
Z+B+1i|-B 


VY Z+«TFL B 
ATRUE FLOOR 
Z+B-11|B 


VY Z+FCL B 
ASCALAR FUNCTION 2+[B 
AUSES:FUZZ FGT TFL FEQ TCL 
Z2+FUZZ<|B 
20 IF~Z 
Z<TCL B 
+0 IF(B FGT TFL B)VB FEQ Z 
Z+TFL B 


VY Z+FPL B 
ASCALAR FUNCTION Z+LB 
AUSES:FUZZ TFL FLT TCL FEQ 
Z+FUZZ<|B 
+0 IF~Z 
Z+TFL B 
+0 IF(B FLT TCL B)VB FEQ Z 
Z+TCL B 


V Z+A FLT B 
ASCALAR FUNCTION Z+A<B 
AUSES:RFUZZ 
Z+A<B-A RFUZZ B 


V Z+A FLE B 
ASCALAR FUNCTION Z+A<sB 
AUSES:RFUZZ 
Z+ASB+tA RFUZZ B 
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V 2+A FEQ B 
ASCALAR FUNCTION Z2+A=8 
AUSES:RFUZZ 
Z+(A RFUZZ B)2\A-B 


Z+A PFGE B 

ASCALAR FUNCTION Z+A2B 
AUSES:RFUZZ 

Z+A2B-A RFUZZ B 


Z+A FGT B 

ASCALAR FUNCTION Z+A>B 
AUSES:RFUZZ 

Z+A>B+A RFUZZ B 


V Z*A FNE B 
(33 ASCALAR FUNCTION Z+A#B 
[2] AUSES: RFUZZ 
[3] Z+(A RFUZZ B)<|A-8 


ADVERB FUNCTIONS 


Scalar functions can be extended to arrays in four ways: Simple 
Scalar Extension, Reduction, Outer Products, and Inner Products. 


Scalar Extension. This function extends scalar functions to arrays 
element by element. The scalar function is represented by the 
function * which may be any of the scalar dyadic functions. SD 
performs conformability and general domain checks, but assumes that 
f will detect most domain errors. 


Vi 2A SD B;CZ;RVLZ;RVLASRVLIB3I 


1 AGLOBAL VARIABLES:£ 
[2] AUSES:F 
C3] "DOMAIN' ERROR(O#TYPE A)aA~Fe'e#=' 
C4] "DOMAIN' ERROR(O#T7YPE B)a~Fe'z='! 
[5] +SINGULAR IF(1=x/pA)V1=x/pB 
C6] "RANK ERROR(p0A)#00B 
C74 "LENGTH' ERRORV/(pA)#p0B 
[8] +L1 IF~0OepA 
[9] 2+A 
[10] +0 
[11] SINGULAR:+(ASINGULAR ,BSINGULAR) IF((1#x/pB)V(ppA)<ppB), 
(1#x/pA)V(ppA)>ppB 
C12} +51 
C13] ASINGULAR:+L2 IF 0#x/pB 
[14] Z+B 
[15] +0 
[16] £22:A+(pB)pA 
C1493, “S64 
[18] BSINGULAR:+L3 IF O#x/pA 
[19] Z+A 
[20] +0 


| ee ae SSMS eS 


[21] L3:B+(pA)pB 

[22] 21:C2Z«x/pA 

(23] RVLZ<«CZp0 

[24] RVLA+,A 

[25] RVLB<,B 

[26] I<0 

[27] LOOP:I<+I+1 

[28] +END IF I>CZ 

[29] RVLZCIJ<RVLALI] F RVLBLT] 

[30] +LOOP 

[31] END:Z+«(pA)pRVLZ 
v 


Reduction. SEDUCTION performs vector reduction. # extends vector 
reduction to higher arrays and performs error checking. The scalar 
dyadic function which is reducing is represented by ie F contains 


the character symbol which denotes the reducing function, and F1 is 
the value of the corresponding identity element. 


F/(I]X is represented here as R I AXIS X with appropriate 
specifications of F, FE, and £1. 


V 2+R B;LZ2;CB;XLB;RVLZ;RVLB;TCB3V3E3Msd3KiI 
C1] AGLOBAL VARIABLES: F £1 
[2] AUSES:OKINDEX SREDUCTION 
[3] I+OKINDEX B 


C4] 'DOMAIN' ERROR(O#TYPE B)a~Ec'#=' 
[5] +L1 IF 0#ppB 

[6] B+,B 

{7] I+,1 


C8) B1:L2+«(( 1+72)+p8),1+p8 
[9] XLB+(p8)(TI 

[10] 7+L2 IF 0#XLB 

[11] 'DOMAIN' ERROR Fe'evno! 
[12] Z+LZpF1 4 
[13] 7+L3 

C14] L2:2<LZpB 

[15] L3:CB<x/pB 

[16] 70 IF(0=CB)v1=XLB 
C17] RVDZ+.Z 

[18] RVLB+,B 

[19] 2CB+«x/I+pB 

[20] V<TCBx 1+i\XLB 

[21] E+TCBx 1+XLB 

[22] M+J+1 

[23] OUTER:K<1 

[24] INWER:RVLZ(M]+SREDUCTION RVLBLJ+V] 
(25] K+K+1 

[26] Jd+J+1 

[27] M<+M+1 

[28] -+IWNER IF KsTCB 

[29] d+Jt+E 

[30] ~+OUTER IF JsCB 

[31] Z+LZpRVLZ 


ADVERB FUNCTIONS 19 


V Z+SREDUCTION B3d 
cag AUSES:F 
[2] J+(10)ppB 
C3] Z+B( J] 
C4] Lid+J-1 
C5] +0 IF J=0 
[6] Z+BL J] FZ 
as +f 


Outer Product. * represents the scalar dyadic function in A °.F B, 
and £ contains the character symbol which denotes the function. 


Vo ZeA OF BS CASCR SC 
AGLOBAL VARIABLE 
AUSES:F 
"DOMAIN" ERROK(O#TYPER A)a~Pe '#=! 
"DOMAIN ERROR(CO#TYP? B)A~Fe te=' 
CA+x/pA 
CB<x/pB 
CZ<CAxCB 
RVLZ+CZp0 
+END IF 0=CZ 

] RVLA+,A 

] RVLB+,8B 

] d+H+0 

] OUTER: d+d+1 

] -+END IF J>CA 

J] K+0 

] INNER: K+K+1 

J 

J 

J 

] 

J 

] 


sRVLZ3RVLASRVLB 3d 33K 


+I NCREMENT IF K>CB 
RVLZUM+KIJ+RVLAL J] F RVLBCK] 
*+INNER 

INCREMENT 3 /4+M+CB 

+OUTER 
END:2+((pA),pB)oRVLZ 


- 
v 


Inner Products, Inner products of the form A F.G B and decode 
expressions 418 are closely related and are combined here. The 
main function, BASEPROD, distinguishes inner product from decode by 
the value of the global variable I, which is set in D#&CODE and IP. 
BASEPROD in turn uses f to reduce or SDECODE for 1 between 
conforming vectors. F and G denote the first and second inner 
product functions, and F is the identity element of F. 


V Z+A IP B 
] AGLOBAL VARIABLES:I 
] AUSES:BASEPROD 
J I+2 
] Z+A BASEPROD B 


VY Z+A DECODE B 


C1] AGLOBAL VARIABLES:I 
[2] AUSES:BASEPROD 

By I+1i 

[4] Z+A BASEPROD B 
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C1] 
(2] 
C3] 
C4] 
C5] 
C6] 
[7] 
C8] 
[9] 


C4 

[2] 

[3] 

C4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 
L134 
[13] 
C14] 
C15] 
{16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
(26] 
{27] 
[28] 
[29] 
[30] 
BK 
(32) 
{33] 
[34] 
[35] 
[36] 
{37] 


V 


Z<A SDECODE B;U;I 
U+(pA)[pB 

A+UpDA 

B<+UpB 

I+1 

2+BL1] 
L:I*«I+1 

+0 IF U<I 
Z2<BCLI1+2xACI] 

+L 


Z+A BASEPROD B;LLA;FLB3L2Z;CZ;RVLZ;VA;TCB;VB;RVLASRVLB 3d 3MSKs 


AGLOBAL VARIABLES:I Fi G@ 

AUSES:G R SDECODE 

I+I 

I<i0 

+IP1 IF I>1 

‘DOMAIN' ERROR O#TYPE A 

'‘DOMAIN' ERROR O#TYPE B 

+L1 
IP1:'DOMAIN' ERROR(O0#TYPE A)a~Ge'==' 
'DOMAIN' ERROR( O#TYPE B)a~Ge'z=' 
L1:7L2 IF O#ppA 

A+,A 
L2:7*7L3 IF Oz#ppB 

B+,B 
L3:LLA+ 1+pA 

FLB<+1+pB 

'LENGTH' ERROR(12LLA )a(12FLB)ALLA#FLB 
LZ+( 149A4),1+08 
ITP2:RVLZ<(x/LZ)pF1ixI>1 
L4:+END IF 0=x/(pA),pB 

VA+ 1+1LLA 

TCB<+x/1+08 

VB<TCBx 1+. FLB 

CA+x/pA 

RVLA<,A 

RVLB<+,B 

J+~M+0 
OUTER:K+1 
INNER:2IP3 IF I>1 

RVLZUM+K J+RVLALJ+VA] SDECODE RVLBCK+VB] 
>L5 
ITP3:RVLZCM+K1+R RVLALJ+VA] G RVLBL K+VB] 
L5:K+K+1 

*+INNER IF K<TCB 

M+M+TCB 

J<+J+LLA 

+OUTER IF J<CA 


[38] END:Z<LZ2pRVLZ 


v 


r 
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MIXED FUNCTIONS 


Mixed functions generally operate on the structure of their 
arguments rather than on the values of the elements. Some mixed 
functions can be written with a subscript in expressions of the 
form A F[I] B and F[I] 8. We represent these expressions, using an 
auxiliary function, as A F I AXIS B and F I AXIS B respectively. 
Eliding I AXIS is analagous to eliding [J]. 


Catenation. A,[7]B is represented by A COMMA I AXIS B 


VY Z*+A COMMA B;L;3;I 
AGLOBAL VARIABLES:I 
AUSES:CAT COMMACHECK 
I+A COMMACHECK B 
LAMIWHATE:*CATEWATE IF T=tI 
+(ASCALAR,BSCALAR) IF(0=ppA),0=ppB8 
"RANK" ERROR(ppA) #008 
"LENGTH ERRORV/(pA)#peB 
>BSCALAR 


t 
- 


CALAR:L+«pA 
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VY 2+A COMMACHECK B 
AGLOBAL VARIABLES :I 
Z+I 
I<+10 
"'DOMAIN' ERROR(TYPE A)#TYPE B 
"INDEX" ERROR O#TYPE Z 
"INDEX" ERROR 250902 
"ITNDEX' ERROR 1<x/p2Z 
+L2 IF 0#x/pZ 
2*+0.5 
] 70 IF(0=ppA)A0=p0B 
] Z+(ppA)Tpo8 
yp 0 
] 
] 


WHF O WW bt ts Lt 


L2:'INDEX' ERROR Z<0 
'TWDEX' ERROR Z21+(ppA)[ppB 


PSS SSS eA a See 
rPrRrRrRrF WW ONOWFwWhH 


= 


V 


V Z2+A CAT B3R3LZ2;NOTT;LA3LB;RD;WA3WB3CZ3RVLZ;RVLA;RVLB;TCZ3VA; 
VB sd 3K3M3I 
AGLOBAL VARIABLES:I 
I+I 
I+«10 
R+(ppA)fppB 
LZ+Ro0 
NOTI<«((2-1)4.1R),I4+iR 
LA+pA 
LB<pB 


meme eSSe 
OnNOn FwWNHR 
eee ee 
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(9) +(ASCALAR,BSCALAR) IF(0=ppA),0=908 
[10] RD+(ppA)-ppB 

[11] ‘RANK' ERROR 1<|RD " 
[12] +(RD1i,RD1) IF(RD=1),RD= 1 

[13] ‘LENGTH' ERRORV/(LA#LB)CNOTI) 

(14) L2Z2CNOTIJ<LACNOTI) 

[15] WA+LA(TI] 

{16] WB+LBCT] 

(17) 751 

[18] RD1i:'LENGTH' ERRORV/LALNOTI1#LB 
[19] LZCNOTII+LB 

[20] WA<LALTI] 

[21] WB+1 

[22] ~Li 

[23] RD1:'LENGTH' ERRORV/LA#LB(NOTI] 
C24] LZCNOTII<LA 


[25] WA+1 

[26] WB<«LB(T] 

C27) +52 

[28] ASCALAR:LZCNOTI]+LA+LB(NOTI ] 
[29] WA+1i 

[30] WB+LB(I] 

C3n3 bt 


[32] BSCALAR:LZCNOTI]+LB<LACNOTI] 
[33] WA«LA[TI] 

[34] WB«1 

[35] £1:L2CI]+wWA+WB 
[36] 72L2 IF O#WA 

[37] 2Z+LZpB 

[38] +0 

[39] L2:*L3 IF O#wWwB 
[40] 2+LZpA 

C41] +0 

C42] 23:CZ+«x/L2Z 

[43] -7L4 IF 0#CZ 

[44] Z<LZpTYPE A 

[553] +0 

[46] L4:RVLZ+CZpTYPE A 
[47] RVLA+(x/LA)pA 
48] RVLB<(x/LB)pB 
{49] TCZ+x/I+LZ 

[50] WA+WAxTCZ 

[51] WB<WBxTCZ 

[52] VA+ 1+.1WA 

[53] VB+ 1+.WB 

[54] J<«K+M+1 

[55] LOOP: RVLZ(M+VA]+RVLALJ+VA] 
(56) M+<M+WA 

[57] J+d+WA 

[58] RVLZCM+VBI+RVLBCK+VB] 
[59] M+<M+WB 

[60] K+K+WB 

{61] -+LOOP IF M<CZ 
[62] Z<LZpRVLZ 
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Compression. A/[IJB is represented by A COMPRESS I AXIS B 


V Z+A COMPRESS B;XLB;LZ;TCB3V3;CZ;RVLZ;RVLB3d3M3K3I 
[ij AUSES:OKINDEX 
C2] I«OKINDEX B 
[3] +L1 IF OzpoB 
C4] B+(I+,1)pB 
[5] L1:+TEST IF 1#x/pA 
[6] A«(pB)([IIpA 


C7] +12 
[8] TEST:'RANK' ERROR 12#ppA 
[9] "LENGTH" ERROR(pB)(IIJ#pA 


[10] L2:"DOMAIN' ERROR~A/Ae 0 1 
[11] L2Z<«pB 
(12] LZCIJ]++/A 
C13] 7+L3 IF(0#+/A)A( (pA) #+/A)A0#x/pB 
Ci4] Z+LZoB 
[15] +0 
[16] L3:7CB+x/I+pB 
{[17] V+ 14+17CB 
[18] CZ<«x/L2Z 
[19] RVLZ+CZoTYPE B 
[20] XLB<+(pB)CTI 
[21] RVLB+,B 
[22] J+M<+1 
[23] OUTER:K+1i 
[24] IWNER:+SKIP IF~ALK] 
[25] RVLZCM+VI+RVLELJ+V) 
[26] M+M+TCBR 
[27] SKIP:K+K+i 
[28] J<+J+TCB 
[29] +INNER IF K<XLB 
[30] +OUTER IF M<sCZ 
[31] 2Z<«LZpRVLZ 
Vv 


Expansion. A\[{I]B is represented by 4 EXPAND I AXIS B 


VY Z+A EXPAND B;L2;TCB;V;CZ;RVLZ;RVLE LAST SMS K3I 


[1] AUSES: OKINDEX 
[2] I+OKINDEX B 
[3] >+L1 IF O#ppB 
C4] B+(I+,1)0B8 
C5] Li:+?TEST IF 1#x/pA 
[6] A+,A 
(7] +L2 
E [8] TEST:'RANK' ERROR 1#ppA 
[9] L2:"DOMAIN' ERROR~A/Ae 0 1 
[10] ‘LENGTH’ ERROR(pB)(I1#+/A 
[11] L2Z<«pB 
[12] L2CI]<«pA 
J (13] +23 IF(pA)#+/A 
[14] Z+B 
(15] +0 


[16] B3:+L4 IF(0#+/A)A0#x/pB 
[17] Z+LZ2pTYPE B 
[18] +0 
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C19] L4:7CB<+x/T+pB 

[20] V+" 1+17CB 

[21] CZ+x/LZ 

[22] RVLZ<CZpTYPE B 
[23] LA+«pA 

[24] RVLB+,B 

[25] J<+M+1 

[26] OUTER:K+1 

[27] INNER:+SKIP IF~A[K] 
[28] RVLZOM+VJ<RVLBLJ+V) 
[29] J+J+TCR 

[30] SKIP:K<+K+1 

[31] M<M+TCB 

[32] +INNER IF Ks<LA 
[33] +OUTER IP Mscz 
[34] 2«LZoRVLZ 


Deal. A?B is represented by A DEAL B 


Two algorithms are used, The first, lines) 13 through 19, requires 
A iterations and 2 words of storage. The second, lines 20 through 
25%; requires at least A iterations but only A words of storage. 
The decision on line 10 reflects the relative costs of the two 
algorithms in APL\360. 


V 2+A DEAL B;I3d 
Ci] AGLOBAL VARIABLES:0 
[2] AUSES:ROLL 


Usd "RANK' ERROR i12#x/pA 
C4] "RANK' ERROR 12x/pB 
[5] "DOMAIN' ERROR O#TYPE A 
[6] "DOMAIN' ERROR(A<0)VA#LA 
[7] "DOMAIN" ERROR O#TYPE B 
[8] ‘DOMATN' ERROR BzelB 
C9] "DOMAIN' ERROR A>B 


[10] -+SHORT IF A<|B+16 
[11] 2+(0-1)4.1B 
[12] ~+END IF A=0 
[13] I+0 
[14] LOOP:J+1+I+(ROLL B-I)-O 
[15] J<F+1 
C46] ZC1,J]<Z(7,T) 
[17] +ZLOOP IF Aor 
[18] END:2<A+zZ 
[19] +0 
[20] SHORT:Z<«10 
[21] OUTER:+0 IF A=pZ 
[22] INNER:I<ROLL B 
[23] -+INNER IF TIeR 
C209 ‘Yen. iF 
[25] -+OUTER 

Vv 
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a 


Matrix Division. HB and ABB are represented by MMD B and A DMD B 
respectively. 


If B is a non-singular matrix, then Z+A# is such that 4=8+.«Z. If 
B is over-specified, then Z is a least squares solution. 8 is the 
Matrix inverse of B. The function is more completely described in 
"The Solution of Linear Systems of Equations and Linear Least 
Squares Problems in API,", M.A. Jenkins, Philadelphia Scientific 
Center technical report number 320-2989. 


V Z+MMD B 
ap | Z+((11+9B)°.=11+98 BB 
Vv 


V Z2<A DMD B3P3LA2;L823;FP3IT3d3M23;223M13;21;SIGNA ;ALFA;U 
1) '"DOMAIN' ERROR O#TYPE A 
2) 'DOMAIN' ERROR O#7YPE B 
J ‘RANK’ ERROR 22008 
J ‘RANK’ ERROR~(ppA)e 1 2 
] "LENGTH' ERROR(1+pA)#1t+oB 
] "LENGTH' ERROR(1+pB)<14+pB 
] LA2+,1 
ce] +ON IF 1=ppA 
{3] LA2<1i+pA 
[10] ON:LB2<«1+pB 
“{11] +AHEAD IF(0#LA2)A0#LB2 
£12] 2Z+«(LB2,LA2)p0 
[13] +FIN 
(14] AHFAD: P<11+pB 
C15] Peetf/LiJl|BtQ(OoF pf /18 
[16] B+Bx(pB)pF 
[17] B«B,A 
[18] I+«0 
[19] LOOP:J+I 
[20] I+f+1 
[21] +END IF LB2<I 
[22] M2<«f/L1]|(0,-LA2)+(J,7)+B8 
[23] ‘DOMAIN' ERROR FUZZ2f /H#2 
C24] F2«+J+M2i[ /M2 
£25] PCI,I2)<+PCI2,7) 
£964. BEST LAI l 322.7) 
[27] Mi+|J+BC 37] 
[28] Ti+d+M1iif /M1 
(29] BCI,713;)<BCI1,7;) 
[30] SIGMA<+/(J+BL3;I])*2 
[31] ALFA+( 1*0<BLI;I])xSIGMA*0.5 
(32] U«BLI;I1-ALFA 
[33] Bld +1itpBsl+1it+pBl<((J,7)+8)-(U,I+BL3;I])°.«( sSIGNA-BCI;I]~x 
ALFA )x(U,I+BL3I1)+.x(J,T)+8 
[34] BCI;I]<ALPA 
[35] +ZLOOP 
[36] END:2Z<(2LB2,2A2)p0 
£37] I+#(10)p1+2LB2 
[38] QBACK:I<I-1 
[39] +RE IF O0=1 
C40] ZO l3]<((LB2+BC13;1])-(LB2+BCI;)])+.x2)+B(I3T] 
[41] +QBACK 
[42] RE: 2+ZCAP;1*x8(OpZ)pF 
[43] FIN:70 IF 1#ppA 
C44] Z+,Z 
Vv 


[ 
[ 
[3 
C4 
[5 
[6 
{7 
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Take and Drop. AtB and A+B are represented by A 7AKE B and 
A DROP B respectively. 


V Z2+A TAKE B 
J AUSES:TAKECHECK TAKER 
2] TAKECHECK 
J Z+A TAKER B 


V Z+A DROP B 
a a USES: TAKER TAKECHECK 
Pes TAKECHECK 
[3] A+(— 1*0<A)xO[(pB)-|A 
[4] Z+A TAKER B 


V TAKECHECK 


i AGLOBAL VARIABLES:A B 
[2] "DOMAIN" ERROR Oz#TYPE A 
3] "RANK' ERROR 2<90A 

C4] A+,A 

[5] "DOMAIN' ERRORV/Az#LA 
C6] >LZ1 IF O=ppB 

{7} "LENGTH' ERROR(ppB)#pA 
[8] +0 

C9] L1:B+((pA)p1)oB 


Vv 


V 2A TAKER B3L23LB;QB;Q23;L;C;BI;2Z133;SB3;S2;SL3;P:RVLB;RVLZ 
[19 A+,A 
[2] LZ+|A 
(33 LB<+oB 
C4] +L1 IFv/LZ*LB 
C5] Z2+B 
{6] +0 
C7] 2i:2«LZo07YPE B 
[8] +0 IF(0=x/LZ)V0=x/LB 
[9] QB+(A<0)xO[LB-L2Z 
[10] QZ+(A<0)xO[LZ-LB 
[ei -LebZ Le 
(12] C+x/L 
(13] BI+ZI<«Cp1i 
[14] J+ppB 
[15] SB+«SZ<SL+1 
C16] LOOP: P<LCJ]|l(~14+1C)+SL 
C17] BI+BI+SBxQB[J]+P 
[18] 2Z27+«2I+SZxQ2ZC71+P 
[19] SB+SBxLB[J] 
[20] $2+SZxLZ[J] 
(21] SL<«SEIxLE7] 
[22] dJd+J-1 
[23] +LOOP IF J>0 
[24] RVLZ+,Z 
[25] RVLB+,B 
[26] RVLZCZI1«<RVLBCBII 
[27] 2Z+LZpRVLZ 
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Encode. ATE is represented by A ENCODE B.- 


V 


[iJ AUSES:SENCODE 
i 2") 'DOMAIN' ERROR O#TYPE A 
[3] 'DOMAIN' ERROR O#T7YPE B 
C4] LZ+(pA),pB 
C5] CZ+x/LZ 
[6] CB<x/pB 
[7] +L1 IF 0#CZ 
[8] Z+LZp0 
[9] >0 
[10] Li:*L2 IF O#ppA 
11] A+,A 
[12] L2:E+x/1+pA 
[13] VA+Ex i+11tpA 
[14] VZ<+CBxVA 
[15] RVLZ+CZp0 
[16] RVLA+,A 
{(17] RVLB<,B 
[18] J+~M+0 
[19] OUTER:K+1 
[20] INNER: RVLZUM+K+VZI+<RVLALJ+VA) SENCODE RVLB(K] 
E224) KekK+2 
[22] ~+INNER IF KsCB 
[23] M+M+CB 
C24] d<«J+1 
[25]  7+OUTER IF J<E 
[26] Z<«LZpRVLZ 
C27] +0 
Vv 
VY Z<A SENCODE B;I 
{iJ Z+(I+pA)p0 
[2] L:+REM IF ACLI1]=0 
{33 Z(IJ+ACII|B 
C4] +0 IF J=1 
C5] B+(B-Z(1I])+ACLT] 
C6] +0 IF 0O=B 
C7] I+I-1 
(8] +L 
{9} REM:ZLIIJ+B 
v 
Membership and Inverse Index. AeB and A.B are 
A MEMBER B and A XOF B respectively. The function FE 
line 15 of xOF to indicate a fuzzed comparison. 
VY Z+A MEMBER B 
C13 AGLOBAL VARIABLES:0 
{2] Z+((p,B)+0-1)2(,B)14 


v 


Z<«A ENCODE B3L23;CA;CB3E3VA3;VZ;RVLZ;RVLA ;RVLB 3d 3M3K 


represented by 
@ is used on 
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VY Z<A XOF B3LA;CB;RVLZ;RVLB3d 3K 
Es AGLOBAL VARIABLES:0 
eS AUSES:FEQ 
£3) "RANK' ERROR 1#ppA 
C4] +L1 IF(0#pA)A0#x/pB 
[5] Z«(pB)p0 
[6] +0 
C7] £1:LA+pA 
[8] CB+x/oB 
C9] RVLZ<CBeO 
[10] RVLB+,B 
[11] J<0 
[12] OUTER:J+J+1 
[13] 7+END IF J>CB 
[14] K+«0 
[15] IWNER:K+K+1 
16] -+EXTRA IF K>LA 
[17] 7+IWNER IF~ALK] FEQ RVLB(J] 
[18] EXTRA:RVLZLd ]+K+O0-1 
[19] -+OUTER 
[20] END:Z+«(pB)pRVLZ 
Vv 


Index Generation. 18 is represented by XGEN B. 


V Z2+XGEN B3d 
C1] AGLOBAL VARIABLES:Q 


(2] '‘DOMAIN' ERROR O#TYPE B 
C3] '"RANK' ERROR 1#x/pB 

C4] 'DOMAIN' ERROR BeL|B 
[5] Z+Bpo0 

[6] J+0 


[7] L:d+J+1 
[8] +0 IF J>B 
[3] ZLJ]+J+Q-1 
[10] +2 

Vv 


Transpose, ®B and AQB are represented by MTRANSPOSE B 


A TRANSPOSE B respectively. 
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and 


VY Z+MTRANSPOSE B3LB3U3d3L23;W;CZ;73;RVLZ;RVLB;RZ3531 5K 


C1] +L1 IF 1<x/pB 
[23 Z+B 

[3] +0 

C4] L1:LB<pB 

[5] U+i+J+poB 

C6] LZ<«W+JpCZ+1 
EF} LOOP: LZ2(U-dJ]+LB[J] 
[8] T+CZxLBld] 

(9] WLU-d ]+T+CZ 
fi0d -€2<7 

(iid J<J=12 

[12] +LOOP IF J>0 
(13) W<W-CZ 

[14] RVLZ+CZpTYPE B 
[15] RVLB+,B 

[16] RZ<«ppB 

(17] S+RZpI+«d+1 
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[18] M 
[19] 
[20] s 
[21] 
[22] 
{23] 
[24] 
25] B 
[26] 
{27] 
[28] 


[9] 2 


[16] ZL 


£207) DB 


C29] 2£ 


[32] LZ 


Ravel. 


v 
C1] 
Vv 


AINLOOP:K+RZ 
RVLZCdJ+RVLB(TI 
EEK:SCKJ<1+S([K]J 
+BACKUP IF LZ(K1<S(K] 
I+I+WLK] 

J+J+1 

+MAINLOOP 
ACKUP:S[K]+1 
K+K-1 

+SEEK IF 0<K 
Z+LZpRVLZ 


Z+A TRANSPOSE B;:LB3RB3RZ2303W13W23LZ23M3CZ23;CL3I;BI 
AGLOBAL VARIABLES:0 
"DOMAIN' ERROR O#TYPE A 
"RANK' ERROR 2<ppA 

A+,A 

"LENGTH' ERROR(ppB)=pA 
+*L1 IFv/AzippB 

Z+B 

>0 

1:'DOMAIN"' ERRORV/(A#LA)VA<O 
A+A+1-0 

RZ<[T/A 

"DOMAIN' ERRORV/~(18Z)eA 
LB+poB 

RB<+p08 

W1i<«(J+RB)pi 
OOP1:WilJ-1])«W1i[J]xLB[J] 
J+J-1 

+L00P1 IF i<d 
LZ+W2+RZp~J+1 

OOP2:M+J=A 

L2CIIJ+l /M/LB 
W20LI1++/M/W1 

J+J+1 

*+LOOP2 IF JsRZ 

CZ+x/LZ 

+~£2 IF 0#C2Z 

Z+LZoTYPE B 

+0 

2:CL+ 14+1CZ 

I+RZ 

BI+CZp0 

OOP: BI+BI+W2CIIxLZ(I1|cL 
CL<LCL:LZ(T] 

I+I-i 

+LOOP IF 0<I 
Z2+LZo(,8)C1+BI] 


»B is represented by RAVEL B. 


Z+RAVEL B 
Z+(x/pB)pB 
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Rotate. A¢([I]B is represented by A ROTATE I AXIS B. 


V Z+A ROTATE B;CB;XLB;NXLB;TCB;V;E;RVLZ;RVLB;RVLA 3d 3M;K3I 
[1] AUSES:SROTATE OKINDEX 
[2] I+OKINDEX B 
Boi 'DOMAIN' ERROR O#TYPE A 
[4] "DOMAIN' ERROR 1€A#LA f 
[5] CB+x/pB 
C6) XLB+(pB)(I]J 
C70 WXLB<«((~1+I)+pB),I+oB 
[8] +TEST IF 1#x/pA 
[93] A+NXLBpA 
(10) +Li1 
(11] TEST:'RANK' ERROR(pNXLB)#ppA 
(12] ‘'ZLENGTH' ERRORV/NXLBxpA 
C13] £21:+*L2 IP(1¢0#XLB|A)aA(1#x/XLB)AO#CB 
[14] Z+B 
[15] +0 
[16] 22:TCB+x/I+pB 
[17] V«+TCBx 1+.XLB 
(18] E«TCBx 1+XLB 
[19] RVLZ<CBpTYPE B 
[20] RVLB+,B 
[21] RVLA+,A 
[22] J+M<+1 
[23] OUTER:K<+1 
[24] INNER: RVLZCJ+VJ+RVLALM] SROTATE RVLBLJ+V] 
[25] K+K+1 
[26] dJ+J+1 
[27] M+M+1 
{28] +IWNER IF KsTCB 
[29] J+Jd+E 
[30] -+OUTER IF J<SCB 
[31] 2+«(pB)pRVLZ 


V Z+A SROTATE B;D;M 
[1] D+pB 
[2] M+D\A 
E3] Z+BL (M4+1D-M),i1M] 


Index check. The following function is used by mixed indexable 
functions to check the validity of an index and to supply an 
assumed value if the index is elided. 


V Z+OKINDEX B 


[1] AGLOBAL VARITABLES:I 4 
[2] <LI 

C3] I+i0 

C4J "INDEX' ERROR O#TYPE Z 

£5] "INDEX' ERROR 25p02 

[6] Z<,Z 


[7] +L IF 0<pZ 

Cs] +0 IF O=ppB 

[9] Z2<+ppB 

[10j] +0 

C11] Z2Z:'INDEX' ERROR 1<pZ 

Ci2a "INDEX" ERROR~ZeippB 
Vv 
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MONADIC INDEXABLE MIXED FUNCTIONS 


The monadic functions 4 ¥ $ and scan are all similar in that they 
are indexed and they preserve the structure of the argument. 
MINDEXED is the common control function. 


VY Z+A MINDEXED B;CB;XLB;RVLZ;RVLB;TCB3V;E3d3K3I 
fe AUSES:SREVERSE GRADE SSCAN OKINDEX 
[2] I*+OKINDEX B 
[3] Z2<B 
[4] CB<+x/pB 
[5] +0 IF O=CB 
[6] +f1 IF Ae 1 4 
E79 Z+(pB)o1l 
C8] L1i:XLB<(pB){T] 

[9] AFOR SCALAR B XLB=10 

[10] +0 IF 1=*/XLB 

[11] RVLZ+RVLB+,B 

[12] fCB+x/I+pB 

[13] V+TCBx 1+1XLB 

[14] E<+TCBx 1+XLB 

[156J] 2 

[16] OUTER:K+1 

[17] INNER :+( REVERSE ,GRADEUP,GRADEDOWN , SCAN ) IF A=14% 
[18] REVERSE:RVLZ(J+*V]+SREVERSE RVLBLJ+V] 
(19] +L2 

[20] GRADEUP:RVLZ{J+VJ+1 GRADE RVLBLJ+V] 
[21] +bL2 

[22] GRADEDOWN:RVLZ(J+V]*2 GRADE RVLBLJ+V J 
£233 352 

[24] SCAN:RVLZ(J+VJ<SSCAN RVLBLJ+V] 

[25] £22:K+K+i 

[26] dJ+J+1 

[27] 7+INNER IF KsTCB 

[28] d+Jt+E 

[29] +OUTER IF JsCB 

[30] 2Z<(pB)pRVLZ 


Gradeup and Gradedown. ACI1]B and Y{IJ]2 are represented by 
GRADEUP I AXIS B and GRADEDOWN I AXIS B respectively. GRADE 1s 
called by MINDEXED to perform the appropriate vector grade. 


VY Z+GRADEUP B 
[1] AUSES:MINDEXED 
[2] 'DOMAIN' ERROR O#TYPE B 
[3] Z+2 MINDEXED B 


Y Z<+GRADEDOWN B 
[1] aAUSES:MINDEXED 
[2] ‘NOMAIN' ERROR O#TYPE B 
[3] Z+3 MINDEXED B 
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V Z+A GRADE B;U3d;K;C 
C1] a A=1 FOR 4, A=2 FOR ¥ 
[2] AGLOBAL VARIABLES:0 
[3] U+pB 
[4] Z+J+1 
[5] OUTER:d+J+1 
C6] +END IF U<d 
[7] C+BLJ] 
[8] K+d 
[9] INNER:K+K-1 
[10] -+EX IF 0=K 
[11] +INNER IF((BCK]>C),BCLK1<C)(A] 
[12] EX:2+(K+Z),J0,K+Z 
(13) C+BLK+1] 
[14] BCK+1]+BL(J] 
(15] BC7]+c 
[16] -+OUTER 
[17] END:2+(Q-1)+Z 
Vv 


Reversal. o{IT1B is represented by REVERSE I AXIS B. 
called by MINDEXED to perform a vector reversal. 


V Z+REVERSE B 
C2] AUSES :MINDEXED 
[2] Z+1 MINDEXED B 
v 


VY Z+SREVERSE B 
C1] Z+BL (1+pB)-1pB] 
Vv 


Scan. F\{I1B is represented by SCAW I AXIS B. F, 


SSCAN, represents a dyadic scalar function. 


SREVERSE is 


called from 


contains the 


character symbol which denotes F. SSCAN is called by MINDEXED to 


perform a vector scan. 


¥V Z*+SCAN B 
eg AGLOBAL VARIABLES:EF 
[2] AUSES :MINDEXED 
[3] 'DOMAIN' ERROR(O#TYPE B)A~Ec'#=! 
C4J Z+4 MINDEXED B 


¥V Z+SSCAN B3;d;K 
C4] AUSES:F 
[2] J+pB 
C3] Z+B 
[4] OUTER: K+d 
[5] INNER: K+«K-1i 
C6] +L IF K=0 
ky A | Z20J1+Z0K] F Z(dJ] 
C8] +INNER 
C9] L:d+J-1 
[10] -7+OUTER IF J>1 
Vv 
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_———— 


AUXILIARY FUNCTIONS 
AXIS is used to represent function subscripting. 


TYPE returns a space if its argument contains characters, otherwise 
it returns zero. 


IF and ERROR are used for convenience, The third line of ERROR 
(not strictly APL) in APL\360 causes a return to the last level of 
immediate execution. 


F and G are’ used to represent scalar dyadic functions in reduction, 
scan, and inner and outer products. 


V Z+A AXIS B 
C1] AGLOBAL VARIABLES:0 I 


[2] ‘TNDEX' ERROR O#T7YPE A 
[3] I+A+1-0 
C4] 2+B 
v 
V Z2+TYPE B 
[1] AGLOBAL VARIABLES:C 
[C2] Z+' 1 
C3] +0 IFv/,BeC 
C4] 2+0 
v 
VY Z+A IF B 
(1) Z<+B/A 
Vv 


VY A ERROR B 
Gay +0 IF~v/B 
{2] A,'ERROR' 
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INDEX OF FUNCTION DEFINITIONS 


ABS 
ACOSH 
AND 
ARCCOS 
ARCSIN 
ARCTAN 
ASINH 
ATANH 
AXIS 
BASEPROD 


COMMACHECK 
COMPRESS 
COS 
COSH 
DEAL 
DECODE 
DMD 
DROP 
ENCODE 
ERROR 
ETO 
EXP 
EXPAND 
F 

FACT 
FCL 
FEQ 
FFL 
FGE 
FGT 
FLE 
FLT 
FNE 
FUZZ 


G 

GRADE 
GRADEDOWN 
GRADEUP 
IF 

IP 

LN 


LOG 

MAX 

MEMBER 

MIN 
MINDEXED 
MMD 

MPLUS 
MTRANSPOSE 


OKI NDEX 
OP 

OR 
PITIMES 
R 
RAPPROX 
RAVEL 
RECIP 
RES 
REVERSE 
RFUZZ 


SDECODE 
SENCODE 
SHRIEK 
SIGNUM 
SIN 

SINH 
SREDUCTION 
SREVERSE 
SROTATE 
SSCAN 
TAKE 
TAKECHECK 
TAKER 

TAN 

TANH 

TCL 

TFL 
TRANSPOSE 
TYPE 

XGEN 

XOF 


INDEX 
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